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UNIFORM  GEOMETRICAL  THEORY  OF  DIFFRACTION 

by 

p.H,  Pathak 

The  Ohio  State  University  ElectroSclence  Laboratory 
1320  Kinnear  Road 
Columbus,  Ohio  43212 


Keller's  geometrical  theory  of  diffraction  (GTD)  represents  a major  breakthrough  In  solving  a wide 
variety  of  electromagnetic  (EM)  radiation  and  scattering  problems  at  high  frequencies.  In  particular,  the 
GTO  Is  an  extension  of  geometrical  optics  to  Include  a class  of  diffracted  rays  via  a generalization  of 
Fermat's  principle.  These  diffracted  rays  are  initiated,  for  example,  from  geometrical  and  electrical 
discontinuities  In  a scatterer,  or  from  points  of  grazing  Incidence  on  smooth  convex  parts  of  the 
scattering  surface.  However,  being  a purely  ray  optical  theory,  the  original  GTD  fails  wtthfn  the 
transition  regions  adjacent  to  geometric  optical  shadow  boundaries  where  the  diffracted  field  generally 
assumes  Its  largest  value.  This  limitation  of  the  GTO  Is  overcome  via  the  uniform  version  of  the  GTO 
(l.e.,  UTD)  which  requires  the  diffracted  field  to  make  the  total  high  frequency  field  continuous  across 
the  optical  shadow  boundaries.  The  UTO  solutions  for  che  diffraction  by  edges  and  smooth  convex  surfaces 
are  reviewed  In  detail  after  Introducing  the  basic  concepts  of  GTO.  Results  based  on  a few  additional  UTO 
solutions  are  also  presented  together  with  a few  selected  applications  of  these  UTO  solutions  to  predict 
the  EM  radiation  and  scattering  from  complex  structures. 

I.  INTRODUCTION  / 

An  efficient  analysis  of  the  radiation  and  scattering  of  waves  by  objects  which  are  large  In  terms  of 
the  wavelength  can  be  performed  via  high  frequency  techniques.  One  of  the  most  versatile  and  useful  high 
frequency  (HF)  techniques  Is  Keller's  geometrical  theory  of  diffraction  (GTO)  (1,2,31  which  was  developed 
in  the  early  1950s.  The  GTO  constitutes  a significant  extension  of  geometrical  optics  (GO)  In  which  a 
class  of  diffracted  rays  are  Introduced  to  exist  In  addition  to  the  usual  rays  of  GO.  These  diffracted 
rays  are  postulated  via  a generalization  or  Fermat's  principle  with  the  knowledge  that  at  high  frequencies 
diffraction,  like  reflection.  Is  a highly  local  phenomenon,  .lust  as  reflected  rays  originate  from  points 
of  specular  reflection  on  an  Illuminated  surface,  the  diffracted  rays  likewise  originate  from  certain 
localized  parts  on  the  surface;  e.g.,  from  geometrical  and  electrical  discontinuities,  and  from  points  of 
grazing  Incidence  on  a smooth  convex  surface  as  shown  In  Figure  1. 

The  shadow  boundaries  divide  the  space  surroundlnq  an  Illuminated  bony  Into  a lit  region  where  the  GO 
incident,  reflected  and  refracted  rays  are  present,  and  Into  a shadow  region  where  these  GO  rays  are 
absent.  Thus,  the  GO  approach  Is  seriously  In  error  within  the  shadow  region  where  It  predicts  a zero 
field;  this  limitation  of  GO  is  overcome  by  the  GTO  since  the  diffracted  rays  penetrate  Into  the  GO  shadow 
zone  to  entirely  account  for  the  field  therein.  Furthermore,  the  diffracted  rays  can  also  enter  Into  the 
lit  region  and  thereby  provide  an  Improvement  to  GO  In  the  lit  region.  The  total  GTO  field  Is  a 
superposition  of  the  field  of  all  the  GO  Incident,  reflected  and  refracted  rays  together  with  the  field  of 
all  the  diffracted  rays  which  pass  through  the  observation  point.  The  Initial  values  of  the  diffracted 
ray  fields  are  given  In  terms  of  the  diffraction  coefficients  just  as  the  Initial  amplitudes  of  the  GO 
reflected  and  refracted  rays  are  given  in  terms  of  the  reflection  and  transmission  coeff1c(ents. 

r 

Oue  to  the  local  nature  of  diffraction  at  high  frequencies,  the.dl ffractlon  coefficients  can  be  found 
from  the  appropriate  solutions  to  simpler  canonical  problems  which  model  the  geometrical  and  electrical 
properties  In  the  neighborhood  of  the  point  of  diffraction  as  In  the  original  problem.  Consequently,  the 
GTD  provides  an  efficient  high  frequency  solution  to  problems  that  cannot  lie  solved  rigorously.  Thus,  a 
GTD  analysis  of  the  radlaMon/scatterlng  from  complex  shapes  can  be  developed  by  simulating  these 
structures  with  simpler  shapes  that  locally  provide  a sufficiently  accurate  description  of  the  dominant 
reflection  and  diffraction  effects.  The  GTD  can  also  he  useful  in  providing  Information  on  ways  to 
control  the  radiation/scattering  from  different  parts  of  the  structure.  It  Is  Interesting  that  even 
though  GTD  Is  a high  frequency  method.  It  Is  often  found  to  work  for  objects  nearly  as  small  as  a 

wavelength  In  size.  Although  GTD  is  not  a rigorous  method,  It  generally  yields  the  leading  terms  In  the 

asymptotic  high  frequency  solutions  of  diffraction  problems. 

Since  the  GTD  is  a purely  ray  optical  theory.  It  falls  within  the  transition  regions  adjacent  to  the 
GO  shadow  boundaries  where  the  HF  field  generally  undergoes  a rapid  transition  across  the  shadow  boundary 
from  one  ray  optical  form  (n  the  lit  region  to  another  ray  optical  form  In  the  shadow  region. 

Consequently,  the  HF  field  departs  from  a strictly  ray  optical  character  within  the  GO  shadow  boundary 
transition  regions.  This  failure  of  the  original  GTO  can  be  overcome  by  uniform  versions  of  the  GTD  such 

as  the  UTO  [4,5]  and  the  HAT  [61.  In  the  present  development,  the  foco*'.  will  be  on  the  UTD.  Basically, 

the  UTD  remains  valid  within  the  GO  shadow  boundary  transition  regions  where  the  ordinary  GTD  falls,  and 
secondly.  It  reduces  to  the  GTO  outside  these  transition  regions  where  the  latter  Is  Indeed  valid. 

The  GTO  ana  its  uniform  versions  (UTD;UAT)  fail  within  the  regions  of  GO  and  diffracted  ray  caustics. 
Ray  caustics  or  focii  occur  whenever  a family  of  rays  (i.e,,  ray  cong  uencesl  merge  or  intersect  to  form  a 
focal  surface,  or  a focal  line  or  a focal  point.  The  field  near  diffracted  ray  caustics  can  be  described 
with  the  help  of  the  equivalent  current  method  (ECM)  [7,8,9]  in  which  the  GTD  indirectly  provides  the 
strengths  of  these  equivalent  currents  that  radiate  fields  at  and  near  the  caustics.  Away  from  che 
caustics,  the  ECM  usually  reduces  to  the  GTD.  The  ECM  can  in  general  be  used  provided  the  GO  shadow 
boundaries  and  caustics  do  not  overlap.  In  the  latter  situation,  ECM  could  in  some  cases  still  be  used 
but  only  after  significant  modification;  alternatively,  the  physical  theory  of  diffraction  (PTD)  can  be 
employed.  The  PTD  was  introduced  by  Ofimtsev  [101  in  the  Soviet  Union  at  about  the  same  time  as  Keller's 
GTD  was  introduced  in  the  U.S.  The  PTD  requires  an  integration  of  the  asymptotic  HF  currents  on  the 
radiating/scattering  body,  if  the  PTO  integrals  can  be  evaluated  asymptotically  outside  the  confluence  of 
GO  shadow  boundary  and  caustic  regions,  then  it  generally  reduces  to  the  GTD.  However,  in  some  special 
instances,  the  GTO  can  be  made  to  work  without  resorting  to  ECM  or  PTJ  despite  a presence  of  a confluence 
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of  caustic  arid  GO  shadow  boundary  transition  regions.  Away  from  the  special  regions  where  It  may  be 
necessary  to  use  ECM  or  PTO,  It  Is  natural  to  employ  the  more  efficient  GTO/UTD  which  unlike  the  ECM  and 
PTO  requires  no  integration  [11]. 

These  notes  will  deal  mostly  with  the  diffraction  by  perfectly-conducting  surfaces  In  free  space. 

The  GTO  formulation  is  presented  after  briefly  Introducing  the  concept  of  wavefronts,  rays 
and  GO  in  Section  11.  Next,  the  UTD  Is  discussed  and  UTD  expressions  are  given  for  the  two  main 
diffraction  mechanisms;  namely,  for  edge  diffraction  and  diffraction  at  a smooth  convex  surface.  Other 
UTO  solutions  are  not  Included  due  to  space  limitations.  Finally,  a few  examples  lllsutrating  the  utility 
of  UTD  to  analyze  radiation  and  scattering  problems  are  given  In  Section  III.  An  ej“t  time  dependence  Is 
assumed  and  suppressed  In  the  following  development. 

II.  THE  6TD  AND  ITS  UNIFORM  VERSION  ~ THE  UTD 

The  basic  Ideas  of  wavefronts,  rays  and  GO  are  briefly  reviewed  at  first.  Diffracted  rays  which 
exist  In  addition  to  the  GO  rays  are  discussed  subsequently. 

A.  Wavefronts  and  Rays 

A wavefront  Is  an  equlphase  surface.  The  connection  between  wavefronts  and  rays  can  be  made  In 
several  ways.  One  such  procedure  which  Is  based  on  the  method  of  stationary  phase  Is  ^escribed  below. 

Let  rjr')  and  TT(r * ) refer  to  the  electric  and  magnetic  field  Intensities  at  any  point  r'  on  an  equlphase 
(or  wavefront)  surface  S.  The  electric  field  lT(r)  at  a point  P ahead  of  the  wavefront  Is  provided  by  the 
equivalence  theorem  as; 

JkZ  e"^R 

E(r)  - If  ds‘  (-$~)[R  x R x 5s(r*)  + Y0R  x Ms(?‘)]  ~J~  (1) 

In  which  the  equivalent  electric  and  magnetic  surface  current  sources  7S  and  Ms,  respectively  on  S are 


Js(r‘)  * n'  x H(r')  ; Ms(r‘)  « E(r')  x n'  . (2a, 2b) 

The  quantity  ZQ  denotes  the  Impedance  of  free  space,  and  'r0"(*0)"1.  Also,  k represents  the  wave  nunber  of 
free  space.  The  vector  R and  the  unit  normal  vector  n'  to  the  surface  S at  r'  are  shown  In  Figure  2. 

Consider  a rectangular  coordinate  system  chosen  for  convenience  so  that  the  x and  y axes  are  tangent 

to  the  wavefront  at  0.  and  OP  * z|0P|  as  In  Figure  2.  It  Is  noted  thut  n'*z  at  0.  It  Is  generally  true 

that  there  Is  at  least  one  point  0 on  S so  that  OP  ■ n ' | OP | ; however,  for  the  present  development  it  is 
assumed  that  there  Is  only  one  such  point  0.  If  there  are  more  points  on  S with  the  above  property  such 

that  the  n‘  directions  from  those  points  Intersect  at  P,  then  P Is  said  to  be  a focal  or  caustic  point. 

From  the  principle  of  stationary  phase  as  described  for  example  by  Silver  [12],  the  e’^k,!  within  the 
Integrand  of  (!)  oscillates  rapidly  for  large  k to  produce  a cancellation  (destructive  interference) 
between  each  of  the  spherical  wave  contributions  to  P which  arise  from  the  different  elemental  sources  on 

-IkR 

ds1  over  S that  do  not  lie  In  the  Immediate  neighborhood  of  0;  whereas,  e changes  slowly  for  the 
spherical  wave  contributions  to  P arising  from  the  elemental  sources  on  ds'  that  are  In  the  Immediate 
neighborhood  of  0 and  thereby  provide  a constructive  Interference  to  P.  Thus,  at  high  frequencies  (or 
large  k),  the  dominant  field  contribution  to  P comes  from  0 on  S;  this  point  0 Is  sailed  the  "stationary 
point."  Without  details  (which  can  be  found  In  [121),  the  stationary  phase  evaluation  of  (1)  yields  the 
following  contribution  from  the  stationary  point: 


o'')ks  ! I np 


|0P|*s  . 


The  expression  in  (3)  describes  the  continuation  of  the  field  at  0 to  the  field  at  P along  the  highly 
localized  or  “ray"  path  W;  the  field  T(P)  In  (3)  Is  thus  referred  to  as  a ray  optical  field.  Figure  3 
shows  a ray  tube  Interpretation  of  the  energy  transport  along  the  central  ray  Up  as  Indicated  by  (3).  The 
PI  and  pj>  In  (3)  refer  to  the  principal  wavefront  radii  of  curvatures  at  0.  From  Figure  3 one  notes  that 


PI  P2 

(pj+s)(p2+s) 


which  Is  automatically  Implied  In  (3).  The  field  E(P)  at  P has  the  same  polarization  as  the  field  E(0)  at 
0 because  the  ray  path  Is  straight  In  a homogeneous  medium.  The  field  Intensity  In  (3)  becomes  singular 
when  s»-|pj|  or  s-~ I P2 I ; these  points  on  the  ray  path  are  marked  (3-4)  and  (1-2)  In  Figure  3,  and  they  are 
referred  to  as  ray  caustics.  The  actual  field  Is  not  stngulai  at  the  caustics;  clearly  the  simple 
expression  In  (3)  is  therefore  not  valid  at  and  near  the  caustic,  even  though  it  Is  asymptotically 
accurate  away  from  the  caustics.  The  distances  pj  and  s>%  are  also  referred  to  as  caustic  distances.  The 
distance  s is  measured  positive  In  the  direction  of  ray  propagation.  The  caustic  distances  p;  and  p^  are 


positive  If  the  caustics  occur  before  the  reference  point  0 as  one  propagate:;  along  the  ray;  otherwise, 
they  are  negative.  If  pi  and  p?  are  positive,  the  wavefront  Is  convex;  It  "hey  arc  negative,  the 
wavefront  Is  concave.  If  one  of  the  radii  (pj  or  P2)  Is  positive  while  the  otier  Is  negative  then  the 
wavefront  Is  saddleshaped.  If  pj  and  P2  are' negative,  and  if:-s>-|pil,  or  V'-jpjl,  then  a caustic  Is 
crossed  at  (3-4)  or  (1-2)  In  Figure  3,  respectively  so  that  ((pi)/(pi+s)1  r,r  ({p2)/(pp+s))  changes  sign 
within  the  square  root  of  (3).  The  positive  nranch  of  the  square  root  Is  chosen  In  (3)  so  that 


f.  I / p 

P + S°  vp  + s 


. « (rh)  * 0 


and  p*pj  or  p2.  Thus,  a phase  jump  of  «/2  occurs  at  each  caustlng  crossing. 


(5) 


The  field  In  (3)  Is  sometimes  also  referred  to  as  an  "arbltre-v1  ay  optical  field  since  pi  and  p2 
can  be  "arbitrary."  The  geometry  In  Figure  3 Is  referred  to  as  :-i  astigmatic  ray  tube  or  a quadratic  ray 
pencil  because  of  the  quadratic  wavefront  surface  approximation  u 0 chat  Is  used  1r,  the  stationary  phase 
approach  leading  to  (3).  It  Is  noted  that  If  pj  and  P2  become  in-  nice,  then  the  field  In  (3)  Is  that  of 
a plane  wave.  If  a\  or  P2  become  Infinite  then  (3)  Is  a cyllndr ' ;-il  wave  field.  Also,  If  p;=p2  (“finite 
value),  then  (3)  Is  a spherical  wave  field.  Thus,  plane,  cyllndi  -‘I,  spherical  and  even  conical  wave 
fields  are  special  cases  of  an  arbitrary  ray  optical  field;  cler'  ,,  It  follows  that  each  of  these  fields 
Is  also  ray  optical. 


Since  the  wavefront  surface  S In  Figure  3 can  be  assoclat  ■ with  either  an  Incident,  reflected  or 
diffracted  wave,  the  field  expression  In  (3)  therefore  applle'  <gually  to  Incident,  reflected  or 
diffracted  rays.  The  field  Is  polarized  transverse  to  the  ra j (id  the  wavefront  at  P Is  “locally"  plane 
If  ks  Is  sufficiently  large  (as  is  assumed  to  be  true  In  the  '.tlonary  phase  evaluation  leading  to  (3)); 
also,  the  local  plane  wave  relation  between  IT  and  H holds,  n . y: 


H(P)  - Y0  s x E(P)  (6) 

or 

E(P)  • - Zfls  x H(P)  (7) 

In  which  s » OP/ | OP | Is  the  ray  direction. 

B.  The  GO  Field 

The  GO  field  Is  a ray  optical  field.  The  Incident  GO  field  Is  associated  with  rays  directly  radiated 
from  the  source  to  the  field  point,  w en  such  an  loci  i*nf  ray  congruence  strikes  an  object,  It  Is 
transformed  Into  a reflected  ray  congruence.  Since  th*  p esent  notes  deal  mostly  with  scattering  by 
Impenetrable  objects,  there  are  no  transmitted  or  refr.f.el  rays  produced  in  this  case.  The  Incident  and 
reflected  GO  rays  satisfy  Fermat's  principle  which  makes  the  Incident  and  reflected  ray  paths  an  extremal. 
Consider  a plane  wave  Incident  on  a perfectly-conducting  wedge  or  a smooth  convex  surface  as  shown  In 
Figures  4(a)  or  4(b).  The  Incident  rays  are  partly  blocked  by  these  surfaces  creating  the  so-called 
shadow  zone  where  the  Incident  ray  optical  field  vanishes.  The  Incident  shadow  boundary  ISB  In  Figure 
4(a)  and  the  surface  shadow  boundary  SSB  In  Figure  4(b)  divide  the  region  of  space  surrounding  the  wedge 
and  the  convex  surface  Into  a lit  zone  and  a shade''  zone. 

It  Is  important  to  note  that  unlike  the  ennv  nttonal  incident  field  which  Is  defined  to  exist  In  the 
absence  of  any  scattering  objects,  the  GO  Inclde,.'.  ray  field  exists  in  the  presence  of  any  objects  that  It 
might  Illuminate,  It  Is  for  this  reason  that  the  GO  Incident  field  becomes  discontinuous  across  the 

shadow  boundaries  ISB  and  SSB  In  Figures  4(a)  an.'.  4(b).  On  the  other  hand,  the  conventional  Incident 

field  would  not  be  discontinuous  anywhere  outside  the  source  region  which  produced  that  field. 

Henceforth,  the  GO  Incident  electric  and  magnetic  fields  will  be  denoted  by  P and  TP , respectlvel 

The  field  of  the  GO  reflected  rays  that  are  produced  by  the  Illuminated  wedge  In  Figure  4(a)  Is  also 
discont  nuous.  In  particular,  the  reflection  shadow  boundary  (RSB)  delineates  the  regions  of  existence 
and  shadow  for  the  reflected  rays  In  Figure  4(a);  whereas,  the  Incident  and  reflection  shadow  boundaries 
ISB  and  RSB  merge  Into  the  SSB  for  the  convex  surface  in  Figure  4(b). 

Consider  a general  problem  of  reflection  where  an  arbitrary  GO  Incident  ray  optical  field  Illuminates 
a smooth,  perfectly-conducting  curved  surface.  The  astigmatic  Incident  ray  tube  associated  with  the 

Incident  ray  In  the  direction  s'  Is  shown  In  Figure  5.  This  Incident  ray  strikes  the  surface  at  0R  to 

produce  a reflected  ray  In  the  direction  sr.  The  astigmatic  reflected  ray  tube  associated  with  the 

reflected  ray  from  On  Is  also  shown  in  Figure  5.  The  field  Er(P)  at  P which  Is  reflected  from  0„  can  be 

written  via  (3)  as: 


/ pr  pr  r 

?>)  « / : L_ i e-Jks  . (8) 

V {»[  * sr)  hrz  + Sr) 

T -r 

It  is  noted  that  E (P)  In  <8)  Is  given  in  terms  of  E (Or ) at  the  point  of  reflection  Or.  Thus,  the 
reference  point  0 in  (3)  corresponds  to  the  point  Or  in  (8),  The  caustic  distances  p^  and  p2r  associated 
with  the  reflected  wavefront  are  shown  In  Figure  5 along  with  the  reflected  ray  distance  sr  from  Or  to  P. 
The  value  of  ^(Or)  is  related  to  the  incident  field  ImOr)  via  tne  boundary  condition 


! 

|# 


f!  x (r1^)  -+f(oR)) » o . ''  (9) 

Here,  n Is  the  unit'  normal  vector  to  the  surface  at  QR.  It  follows  from  (9)  that 

<(%)  -if  - rt(on>  . <10> 

where  W Is  the  dyadic  reflection  coefficient  of  the  surface  at  Or.  Incorporating  (9)  Into  (8)  yields 

Er(P)  « Fttn)  • R /_! L_  . _!i_  e‘Jl<sr  . (11) 

V pj+sr  p£+sr 

-P 

The  reflected  magnetic  field  H (P)  Is  found  easily  from  (11)  via 

fT (P)  ~ Y0sr  X ^(P)  . (12) 


i _ _ * ApA 

It  Is  convenient  to  express  E (QR)  and  i (P)  In  terms  of  the  unit  vectors  (e,^)  anfj  (e,,eJL)  Wh1  ch  are 
fixed  In  the  incident  and  reflected  rays,  respectively,  as  shown  In  Figure  6.  The  ej.e^  and  s'  are 

ApAAp  *|  *p 

mutually  orthogonal;  likewise,  and  s are  also  a mutually  orthogonal  set.  Furthermore,  e,  and  e,  11 

A 1 A A p 

In  the  plane  of  Incidence  defined  by  s and  n at  QR.  As  a result  of  Fermat's  principle,  s also  lies  In 
1 r 

the  plane  of  Incidence  and  6 "6  In  Figure  6.  Thus,  If 

E1 (Qr)  • e](Qr)  e]  + e|(Qr)  ex  (13) 

and 

Er(P)  * Ef(P)ef  + Ej(P)ex  (14) 

» *1  n aj  n ■ 

where  e^s  ’ x e,’  , then  R In  (10)  subject  to  the  boundary  condition  (9)  becomes 

r - ;;  ej  Rh  ♦ a,  a,  rs  m 

h 

* 

In  matrix  notation,  the  above  R can  be  written  as 


Rh  0 

1 0 

o 1S 

0 -1 

Therefore,  in  matrix  notation,  (11)  becomes 

I I I 1 oil  E|«d  I 


e£(P)  o -1  e{(0r)  V (pj+sr)<P2+sr) 

The  caustic  distances  or  the  principal  radii  of  curvature  of  the  Incident  and  reflected  wavefronts  which 

j j p p 

are  denoted  b-t  (p. ,p-,)  and  (p,,p,),  as  well  as  their  principal  wavefront  directions 
are  given  In  (.131:  1 

It  Is  clear  that  the  00  representation  of  (11)  falls  at  caustics  which  are  the  Intersection  of  the 
paraxial  rays  (associated  with  the  ray  tube  or  pencil)  at  the  lines  1-2  and  3-4  as  shown  In  Figure  3. 

Ip  ip 

Upon  crossing  a caustic  In  the  direction  of  propagation,  (p  * +s  • ) changes  sign  under  the  rad’cal  In 

and  a phase  jump  of  *«/2  results  as  explained  earlier.  Furthern»re,  the  reflected  field  Er  of  (11)  fails 
In  the  transition  region  adjacent  to  SS8  of  Figure  4(b).  It  is  Important  to  note  that  near  the  5SB  (l.e., 

as  pj  and  p£  approach  the  following  limiting  values: 


fj 


r og(0R)  “sS1 


for  81  * v/2 


where  pg(QR)  Is  the  surface  radius  of  curvature  In  the  plane  of  Incidence  at  i)R,  and  ph  Is  the  radius  of 

curvature  of  the  Incident  wavefront  In  the  (t,b)  plane  (l.e.,  In  the  plane  tangent  to  the  surface)  at  QR 

for  8^-m/2.  Furthermore,  the  principal  directions  xj  and  x£  of'  the  reflected  wavefront  approach  the 
following  values  for  grazing  Incidence:  1 

5[*b  (at  Qr)  (19) 

for  81  * x/2  , 

x£  * (-sP  » *[)  ♦ n (at  Ojj)  (20) 

- » * « GO 

where  t Is  the  direction  of  grazing  Incidence  at  0R  and  b«t  x n at  0R.  The  total  00  electric  field  E at 

Pl  In  the  lit  region  Is  the  sum  of  the  Incident  and  reflected  ray  optical  fields;  hence, 


i®°(PL)  ~ E1  (PL)  + ^(Qr)-  R 


(Pl+sr)(p2+sr) 


In  summary.  It  Is  noted  that  the  00  incident  and  reflected  fields  are  discontinuous  across  their 
associated  shadow  boundaries  such  as  ISO,  RSB,  and  SSB  In  Figures  4(a)  and  4(b).  The  failure  of  00  to 
account  for  a proper  non-zero  field  within  the  shadow  region  behind  an  Impenetrable  obstacle  can  be 
overcome  through  the  GTD  and  Its  uniform  versions.  Nevertheless,  00  generally  yields  the  dominant 
contribution  to  the  total  high  frequency  fields,  and  It  constitutes  the  leading  term  In  the  0T0  solution. 

The  reflected  GO  field  ?r(PL)  for  'he  two-dimensional  (2-0)  case  can  be  deduced  directly  from  the  3-0 
case  by  allowing  pj  to  approach  Infinity,  'hus,  one  may  let  pR  s prand  pj  ♦ » In  (11)  to  arrive  at  the 
2-D  reflected  GO  field  Er(PL)  as 

Er(PL)  * E1  (0R)  • R j^prr  8'JkSr  - (22) 

In  which  the  Incident  ray  optical  field  (Oft ) I:  now  a cylindrical  wave  at  Or,  and  the  caustic  distance 

Pr  In  (22)  for  the  2-0  case  Is  given  by 


p'  ’ s'  + Pg<V  ’ <23) 

where  e’  has  the  same  meaning  as  before,  and  s’  is  the  radius  of  curvature  of  the  Incident  cylindrical 
wavefront  at  Or.  If  the  cylindrical  wave  Is  produced  by  a 2-0  line  source,  then  s’  In  (22)  can  be  chosen 
to  be  the  distance  from  that  line  source  to  the  point  of  reflection  Or  on  the  2-0  boundary.  The  quantity 
Pg(Oft)  In  (23)  denotes  the  radius  of  curvature  of  the  2-0  boundary  at  the  point  of  reflection  Or. 

C.  The  Diffracted  Ray  Fields 

The  diffracted  rays  are  Introduced  in  the  GTO  via  a generalization  of  Fermat's  principle  as  stated 
previously.  Away  from  the  point  of  diffraction,  the  diffracted  rays  behave  according  to  the  laws  of  GO. 
The  Initial  value  of  the  diffracted  ray  field  Is  given  In  terms  of  a diffraction  coefficient.  The 
phenomenon  of  edge  diffraction  will  be  discussed  first,  and  It  will  be  followed  by  a discussion  on  the 
phenomenon  of  diffraction  at  a smooth  convex  surface.  The  latter  phenomenon  Is  more  complicated  than  the 
first. 

(1)  Edge  Diffraction 

When  on  Incident  ray  strikes  an  edge  In  an  otherwise  smooth  surface,  It  produces  diffracted  rays 
which  lie  on  a cone  about  the  tangent  to  the  edge  at  the  point  of  diffraction  such  that  the  angle  8R 
between  the  Incident  ray  and  the  edge  tangent  equals  '.he  half  angle  of  the  diffracted  ray  cone  as  shown  In 
Figure  1(a).  This  cone  of  diffracted  rays  Is  sometimes  referred  to  as  the  'Keller  cone,"  and  it  results 
from  the  generalization  of  Fermat's  principle  to  describe  rays  diffracted  by  an  edge. 


dam  Aixaskit«a«V4Sx{W^  v( 
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Let  an  arbitrary  ray  optical  field  be  incident  on  a perfectly-conducting  curved  wedge  as  shown  In 
Figure  7.  The  resultant  total  Hf  electric  field  1T(P)  at  any  point  P exterior  to  the  wedge  Is  given  by 

E(P)  = EG0(P)  + Ed(P)  (24) 

-GO 

where  the  GO  field  component  E (P)  Is  given  as 

iS0(P)  =■  eV)!),  + Er(P)Ur  . (25) 

c -■(  -p 

The  domains  of  existence  of  the  Incident  and  reflected  ray  optical  fields  E (P)  and  E (P)  are  Indicated 
by  the  step  functions  0)  and  Ur,  respectively,  which  are  defined  as  follows: 


“1 


and 


1,  If  0 < $ < ir<V 

0,  If  < $ < n* 


1.  If  0 < <1  < n-y 
0,  If  «-$'  < 4 < n*  . 


(26) 


(27) 


The  azimuthal  angles  4 and  4'  are  made  by  the  projections  of  the  directions  of  Incidence  and  observation 
on  a plane  perpendicular  to  the  edge  at  the  point  of  diffraction  Oj,  These  angles  are  measured  from  a 
plane  tangent  to  the  “0"  face  of  the  wedge  at  Or  as  shown  In  Figure  8.  The  plane  tangent  to  the  other 
face  of  the  wedge  at  Of  is  denoted  by  “n*;“  It  is  also  shown  In  Figure  8. 


The  Interior  wedge  angle  is  therefore  given  by  (2-n)».  The  expressions  for  the  GO  Incident  and 

reflected  fields  have  been  discussed  previously.  The  diffracted  field  Ed  exists  exterior  to  the  wedge 
(l.e.,  for  0 < 4 < n»).  From  (2),  one  may  write  the  general  field  expression  for  theray  diffracted  In  the 

direction  sd  from  Of  as: 


Ed(P) 


Ed(P„) 


4 4 


<*?  * *o>(»2 


°d> 


(28) 


The  diffracted  ray  tube  corresponding  to  (28)  Is  shown  In  Figure  7.  The  superscript  "d"  on  4 and  sd 
denotos  that  these  quantities  are  associated  with  the  diffracted  ray  field  component.  In  ordefto  relate 

Ed(P)  to  the  Incident  field  at  tha  point  of  edge  diffraction  0E,  one  moves  the  reference  PQ  In  Figure  7 

to  the  point  of  diffraction  0E  on  the  edge  by  letting  pd  + 0 so  that 


UP)  ■ Hm  C Ed(P0)l 
Pd  ♦ 0 


d 

»2 


-jks“ 


I d * «dw  d . d,  e 
(Pj  ♦ s0)(o2  + sa) 


(291 


(30) 


Since  E°(P)  Is  Independent  of  the  reference  point  P^,  the  above  limit  exists  and  It  Is  defined  as 

11m  Ed(P0)  a E1  (0£>  • , 

p{  * 0 

where  8k  > (>,+•, SQ ; k)  Is  Keller's  "dyadic  edge  diffraction  coefficient"  which  indicates  how  the  energy 

Is  distributed  In  the  diffracted  field  as  a function  of  the  angles  4,  4',  and  6 ; Sk  also  depends  on  n and 
the  wavanurber  k.  From  (29)  and  (30),  It  Is  clear  that  0 e 


Ed(P)  - e1  (Qf)  • . 


(31) 


where  11m  p«  = p.  (edge  diffracted  ray  caustic  distance),  and  likewise  11m  s|J  = sa,  as  shown  in  Figure 

pV  C p«*  0 

7.  £ (P)  is  polarized  transverse  to  the  diffracted  ray  direction  s since  the  field  Ea(P)  Is  ray  optical; 
thus,  the  associated  magnetic  field  can  be  expressed  as 


0d(P)  ~ \ sd  x Ed(P) 


(32) 


3 


If  the  Incident  field  E1 (Qc)  exhibits  a rapid  spatial  variation  at  Q-  then  an  additional  term  referred  to 
as  a slope  diffracted  field  must  be  Included  in  (31)  to  describe  the  diffraction  effects  accurately; 
however,  that  slope  diffracted  field  will  not  be  described  here.  An  expression  for  finding  the  diffracted 
ray  caustic  distance  p6  Is  given  later  In  (43b). 

It  Is  convenient  to  express  the  dyadic  edge  diffraction  coefficient  0^  In  terms  of  unit  vectors  fixed 

*(  a 

In  the  Incident  and  diffracted  rays  as  follows.  Let  s and  e define  an  edge  fixed  plane  of  Incidence 

where  e Is  the  edge  tangent  at  0^.  Likewise,  let  s^  and  e define  the  edge  fixed  plane  of  diffraction. 

The  law  of  edge  diffraction  which  defines  the  Keller  cone  Is  s'-e  » s^e.  Let  80‘  and  a be  parallel  to 
the  edge  fixed  planes  of  incidence  and  diffraction,  respectively  as  *n  Figure  8,  and  let0 


.i  .. 

Bo  * s x 4 


(33a; 33b) 


Here,  4 and  4'  point  In  the  direction  of  Increasing  angles  4 and  4“,  respectively.  The  Incident  field 

-4  *4  a!  a 

E (Qj:)  can  then  be  expressed  In  terms  of  the  triad  of  unit  vectors  (s  ,60,4')  fixed  In  the  Incident  ray; 

likewise,  the  edge  diffracted  field  Ed(P)  can  be  expressed  In  terms  of  (s'1, 8 ,4)  fixed  In  the  diffracted 
ray.  Thus, 

* *0  * ♦*=!•  (34a) 

and 

tV)  * »0  E80  + ♦ E4  • (34b) 

Then 

“e  ■ *®o  »oDes  °eh  * (34c) 

k k 

The  r>e$  and  0eh  can  be  found  from  the  asymptotic  solutions  of  appropriate  canonical  wedge  diffraction 
problems;  they  are  given  by: 


v^' 


n/HxksInS, 


• v 4-41-'  * » 

ig  cos  „ - co'  (V")  cos  S'  - cos  ( V") 


It  Is  noted  that  the  Keller  edge  diffraction  coefficient  In  (35)  beones  singular  at  the  Incident  shadow 
boundary  (ISB)  and  the  reflection  shadow  boundary  (RSB)  which  occur  when  4»*+4'  and  4»x-4‘,  respectively. 
Thus,  the  result  In  (31)  together  with  (34c)  and  (35)  Is  not  valid  at  and  near  the  BO  Incident  and 
reflection  shadow  boundaries.  This  deficiency  of  the  GTD  can  be  overcome  via  the  use  of  untform 
geometrical  theory  of  diffraction  (UTO).  According  to  the  UTD  [4,51,  the  total  HF  field  exterior  to  the 

wedge  Is  still  given  by  (24)  as  In  Kel’«.  \ original  GTD;  however,  the  E In  (24)  and  (31)  Is  now  modified 

-k  * 

so  that  0e  of  (31)  Is  replaced  by  the  HTO  edge  diffraction  coefficient  he  so  that: 


eV)  * r’  (0E).  ne(4.  4'.  8„;  k) 


V S^P.+S*) 


The  De  In  (36a)  can  also  be  expressed  as 

6 a -B  B D -4*0.  . 

e Vo  es  * eh 


In  matrix  notation,  (36a)  becomes 


' °e  __  e-jks 
d, 

s (p  +s  ) 
e 


in  which  the  0es  and  0ep  are  [131: 


Dffi  (♦.♦*!»„)  - 


col  (— -y--)  FfkL'a  ($-$' )1  + cot  ( ~ ) FfklVtM'H 


2n/2»ks1nSo 

cot  F[kLrna+{*4V )!)  + cot--^---  FrkLroa‘(+n')1 


(38) 


where  the  asymptotic  large  parameter  kl  (with  the  superscripts  1,  rn,  ro  on  L omitted  for  convenience)  Is 
required  to  be  sufficiently  large  (generally  greater  than  3)  and 

t inn/  - (6) 

a (8)  * 2 cos2  ( j ) . (39a) 

+ 

The  N are  the  Integers  which  most  nearly  satisfy  the  equation: 


2«nN*  - 8 ■ tv 


with 


6 » $ * , 


(39b) 


(39c) 


Note  that  n“2  for  a half  plane  or  a semi-infinite  curved  screen.  Also,  n*3/2  for  an  exterior  right  angled 
wedge,  etc. 

+ - 4 

For  exterior  edge  diffraction  N «0  or  1,  and  N *-l,  0,  or  1.  The  values  of  N at  the  shadow  and 
reflection  boundaries  as  well  as  their  associated  transition  regions  are  given  In  Table  I for  exterior  wedge 
angles  (1  < n < 2): 


TABLE  I 


The  cotangent  Is  singular  when 

value  of  N 
at  the  boundary 

cot  ( 2n  ) 

♦ *♦'-*,  an  ISB 
surface  ♦•■0  Is  shadowed 

N+  « 0 

cot  ( 2n  ) 

♦ * $'+*,  an  ISB 
surface  **n»  Is  shadowed 

N"  » 0 

.»+(*♦♦'). 
cot  ( 2n  ) 

♦ » (2n-l)*-$',  an  RSB 
reflection  from  surface  $*ni 

Nl  * 1 

. 

cot  ( 2n  ) 

$ * an  RSB 

reflection  from  surface  ♦*() 

N'  * 0 

L 

For  a point  source  (or  spherical  wave)  type  illumination,  the  distance  parameter  L1  1 s :• 

L'  * sin2fio  ’ (40) 

in  which  s*  and  sri  are  the  distances  from  the  point  of  edge  diffraction  at  Of  to  the  source  and  observation 
points,  respectively.  Only  for  a straight  wedge  with  planar  faces  that  Is  Illuminated  by  a point  source, 


* Lrn  » Lf 


si  sd  g 
si"  S0 


(41) 


as  in  (40).  For  an  arbitrary  ray  optical  illumination  which  Is  characterized  by  two  distinct  principal 
wavefront  radii  of  curvature,  pj  and  p^,  the  above  must  be  modified  as  shown  below  In  the  general 
expressions  for  lro  and  lrn  pertaining  to  a curved  wedge;  thus, 


sd(pg+sd)p|pJs*nZB0 

P^(p}Hd)(p^Sd) 


at  ISB 


(42a) 
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f sd(pg+sd)pjp5s1n2?0 

_ Pe(pJ+sd)(p2+sd)  at  >5B  * ^ 

Here,  lro  and  Lrn  are  the  values  of  Lr  associated  with  the  ”0"  and  "n"  faces  of  the  wedge,  respectively. 
Furthermore,  p£  Is  given  by: 

1 1 2fn»ne)(sd -n) 

Z . (43a) 


Mso  pe  In  (36)  Is  given  by: 
1 i _ ne-(s^-sd) 


The  unit  vector  n Is  defined  In  Figure  8(b);  whereas,  ng  Is  a unit  vector  normal  to  the  edge  which  Is  directed 

away  from  the  center  of  edge  curvature  at  0F.  The  rtdlus  of  edge  curvature  Is  denoted  by  “a“  In  (43).  p'  Is 
the  radius  of  curvature  of  the  Incident  wavefront  at  0^  which  lies  In  the  edge  fixed  plane  of  Incidence.  eIn 

the  far  zone  when  sd»p!  ,,  sd»pj  ,,  and  sd»p.  then  the  L(  and  Lp  In  (42a)  and  (42b)  simplify  to 
P1P2  s1n2P0 

L « In  which  the  appropriate  superscripts  on  L,  Pj  and  p2  are  omitted  for  convenience.  It  Is 

noted  that  U and  Lr  In  (42a)  and  (42b)  are  calculated  on  the  appropriate  shadow  boundaries.  The  transition 
function,  F which  appears  In  (38)  contains  a Fresnel  Integral;  It  Is  defined  by 

F(x)  - 2J/7  eJX  r dTe"JX2  • (44) 

fk 

A plot  of  the  above  F(x)  Is  Illustrated  In  Figure  13.  In  (44),  /x  * |/x|  If  x>0  and  /x  * -J  |/7|  If  x^.  If 
x<0,  then  F ( x ) I « F*  (|x|)  where  * denotes  the  complex  conjugate.  Exterior  to  the  (JSjj)  transition  regions 
| X<0  RSB 

x becomes  large  and  F(x)*l  so  that  the  uniform  Des  In  (38)  then  reduces  to  Keller's  form  as  It  should;  namely, 

eh 

-k 

De  ♦ De,  outside  the  transition  region.  (45) 

Near  the  (ISB  and  RSB)  boundaries,  the  small  argument  approximation  for  F ( x ) may  be  employed  (since  x ■ o on 
ISB  and  RSB);  namely,  one  can  Incorporate 

F(x)  ♦ fix  e J(T  + x)  (46) 

x*o 


into  (38)  to  arrive  at  the  following  result  for  the  diffracted  field  E at  ISB  or  RSB: 

jd  ",  It  continuous  ~| 

t 5 f 2,  r ’ + (higher  order)  , 

IS8;RSB  2 terms  ( 

(47) 

on  lit  side  of  IS0;RSB 
’ ' on  shadow  side  o'  ISB;RSB  * 

The  above  result  In  (47)  ensures  the  continuity  of  the  .otal  HF  fiela  In  (24)  at  ISB  and  RSB.  The  field 
contribution  arising  from  the  edge  excited  “surface  diffracted  rays"  Is  not  Included  in  (24);  ft  may  be 
important  for  observation  points  close  to  the  surface  shodow  boundaries  (SSB)  associated  with  the  tangent  to 
the  “0“  and  “n“  faces  of  a curved  wedge  at  Of  If  the  “0"  aid  “n“  faces  are  convex  boundaries.  The  result  In 
(36a)  and  (36b)  along  with  (38)  Is  valid  away  from  any  dlffi  acted  ray  caustics  and  away  from  the  nd"e  cau'tic 
at  Oj. 

For  grazing  angles  of  Incidence  on  a wedge  with  planar  faces,  l>t-  ’ o,  and  Dep  must  be  replaced  by  (l/g) 
°eh'  The  reason  for  the  '/?  factor  in  the  latter  case  Is  explained  a:  follows.  The  incident  and  reflected  CO 
fields  tend  to  combine  into  a single  "total  Incident  field"  as  one  app  oaches  grazing  angles  of  incidence; 
consequently,  only  half  of  this  “total  field"  illuminating  the  edge  at  grazing  constitutes  the  incident  00 
field  while  the  other  half  constitutes  the  reflected  00  field.  The  cast  of  grazing  angles  of  incidence  at  an 
edge  in  a curved  surface  cannot  be  handled  as  easily  as  the  case  of  a weige  with  planar  faces.  Presently,  one 

! 2 l*/3 

can  only  treat  angles  of  incMence  that  are  greater  than  |_kpo(0fT  I where  pg  (Or)  Is  the  radius  of 
curvature  of  the  surface  in  the  direction  of  the  Incident  tKe  point  of  >.Tqe  3’ffractlon  Oj. 
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Under  the  above  restrictions,  the  result  In  (38)  for  pec  simplifies  In  the  case  of  a plane  or  curved 
screen  (n=2  case)  to  en 


-e'j? 

Des(*'*‘*6o)  ' 2/Sk  sins 
eh  o 


secl^-jFru^ff-*'))  t sec(-^“)F[kLra(fH’'l 


(48) 


where  a(S)-2  cos2(S/2)  and  L*>r  are  in  (42a;h)  with  the  urderstandlng  that  Lr  Is  evaluated  at  the  RSB 
corresponding  to  the  face  which  Is  Illuminated;  hence  the  superscripts  "o“  and  “n“  In  Lr  are  dropped  for  this 
n«2  case. 

The  edge  diffracted  field  Ed (P)  for  the  2-D  situation  can  he  obtained  from  (36a)  by  allowing  pe  to 
approach  Infinity  and  by  requiring  B0  • */2;  thus,  for  the  2-D  case, 


e-jksrt 

Ed(P)  ■ E1  (0E)  • 0(!(v,4',*/2;  k)  . (49) 

The  Oe  in  (49)  for  the  2-f)  case  Is  available  from  (36a)  and  (36b)  with  50-«/2  (or  s1nS0  * 1).  Also,  L1  for 
the  2-D  case  is  given  by  (41)  with  S0  « s/2;  In  particular, 

. 1 s'  sd 

L * a . (501 


Likewise,  Lr  Is  obtained  from  (42b)  with  a * »/2,  pi  ♦ «,  pi  a pr  (as  In  (22)),  ind  — ; therefore,  In  the 

2-0  case  o i < e 


r d 
p s 

Pr  ♦ sd 


(51) 


Note  that  pi"  In  (51)  Is  the  same  as  the  one  In  (23);  however,  pr  Is  In  general  different  for  the  "0“  and  “n" 

faces  of  the  wedge,  with  ino  and  Lrn  denoting  the  values  of  Lr  for  these  two  different  faces.  While  the 

expression  for  Lr  in  (42b)  Is  fixed  to  Its  value  on  the  RSB  for  convenience,  the  one  In  (51)  can  be  evaluated 
as  a function  of  the  observation  point  with  almost  the  same  ease  as  If  one  had  approximated  the  value  of  Lr  by 
Its  value  at  the  RSR.  The  values  of  Li  and  Lr  for  the  3-0  case  Involve  various  caustic  distances  as  Is 

evident  from  (42a)  and  (42b).  These  distances  are  generally  slowly  varying  within  the  1S8  and  RSB  transition 

regions  and  It  Is  t.  erefore  convenient  to  approximate  L1  and  Lr  throughout  the  transition  regions  by  their 
values  at  the  IS8  and  RSB  as  done  In  (42a)  and  (42b).  Outside  the  respective  transition  regions,  the  F 
functions  containing  L*  and  Lr  approach  unity  anyway  unaffected  by  the  above  approximation. 

It  Is  noted  that  the  comment  below  (47)  In  regard  to  grazing  Incidence  Is  also  valid  for  the  2-0  case. 

9 »|( 

It  Is  further  noted  that  the  essential  difference  between  De  and  De  Is  that  the  former  is  range  dependent 
whereas  the  latter  Is  not-  As  a result,  (36a)  Is  not  ray  optical  within  the  ISB  and  RSB  transition  regions; 

* “k 

exterior  to  these  regions,  D„  ♦ Oe  as  Indicated  before.  Figure  9 Illustrates  the  diffraction  of  a plane  wave 
by  a perfectly-conducting  half-plane.  It  Is  noted  that  the  geometrical  optics  field  Is  discontinuous; 
however,  the  UTO  diffracted  Field  cancels  the  GO  discontinuity  to  yield  a total  UTO  field  which  Is 
continuous. 


(11)  Diffraction  at  a Smooth  Convex  Surface 

The  geometry  for  this  problem  of  the  diffraction  by  a smooth  convex  surface  Is  shown  In  Figure  10.  The 
totai  high  frequency  field  t(P)  for  the  situation  in  Figure  10  can  be  written  as 


E(P) 


eVl)u  + EfPk)U  * Ed(PL) 
eVj)  Cl  - U] 


, If  P ■ PL  In  the  lit  zone. 

, If  p * Ps  In  the  shadow  zone. 


(52) 


-1  -r 

The  Incident  and  reflected  fields  E and  E are  associated  with  the  incident  and  reflected  GO  rays  shown  tn 
Figure  11.  The  step  function  H In  (b2)  Is  defined  below  with  respert  to  the  surface  shadow  boundary  (SS8) 
as: 


1,  In  the  lit  region  which  lies  above  the  SS8 

U « ( 

0,  In  the  shadow  region  which  lies  below  the  SS8,  (53) 


The  surface  diffracted  field  E (Ps)  follows  the  surface  diffracted  ray  path  Into  the  shadow  region,  as  in 
-d 

Figure  11;  whereas,  the  field  E (Pl)  which  is  diffracted  Into  the  lit  region  follows  the  reflected  ray  path 
-r 

(of  E ) in  this  solution.  Therefore,  it  is  convenient  In  this  problem  to  combine  the  GO  reflected  field 

E (PL)U  and  the  diffracted  field  E ( Pj_ } Into  a single  "generalized  reflected  fi  d",  E9r(P^)U  in  the  lit 
region  so  that  (52)  becomes 


1 


'Sr?--'" 
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m 


i1(PL)U  + E9fPL)U  , If  P * PL  in  the  lit  zone, 
-d. 


E If,)  [1-1)1  , If  P * Ps  in  the  shadow  zone. 

V 

The  fields  lfr(PL>  and  ^(Pj)  are  given  symbolically  by 

r’Ki  - iVeUi  /HE ~ . 

(p[+Sr)(pg+Sr) 


(54) 


(55) 


E^Vs)  "•  ^(Oii'^bihg  +<®J,nin2] 


Fs  e 

T^TpT+T^T 


-jks 


(56) 


where  the  points  0„  and  0,,  and  the  distances  sr  and  sd  are  Indicated  In  Figure  10.  The  surface  diffracted 
ray  caustic  distance  p 15  shown  In  Figure  11.  The  quantities  within  brackets  Involving#  andjtt  In  (55)  and 
^ n b 

(56)  may  be  viewed  as  generalized  dyadic  coefficients  for  surface  reflection  and  diffraction,  respectively. 

It  Is  noted  that  (55)  and  (56)  are  expressed  in/arlantly  In  terms  of  the  unit  vectors  fixed  In  the  reflected 

and  surface  diffracted  ray  coordinates.  The  unit  vectors  ej,  e^,  and  eA  In  (55)  have  been  defined  earlier  In 

connection  with  the  reflected  field.  It  can  be  shown  that  cross  terms  actually  exist  In  the  above  generalized 
dyadic  reflection  coefficient;  hut,  In  general  their  effect  Is  seen  to  he  weak  within  the  SSB  transition 
region.  Also  these  terms  vanish  In  the  deep  lit  region  and  on  the  SSB,  hence  they  have  been  9nored  In  (55). 

At  Oj,  let  tj  be  the  unit  vector  In  the  direction  of  Incidence,  nj  be  the  unit  outward  normal  vector  to 

A * * A A A 

the  surface,  and  bj  « tj  xnj ; likewise  at  Q2,  let  a similar  set  of  unit  vectors  (tg.  n2,  h?)  be  teflned  with 

t2  In  the  direction  of  the  diffracted  ray  as  In  Figure  12.  In  the  case  of  surface  rays  with  zero  torsion, 

bj  * b2.  It  Is  clear  from  Figure  11  that  p$  In  (56)  Is  the  wave-front  radius  of  curvature  of  the  surface 

diffracted  ray  evaluated  In  the  b.  direction  at  02.  First,  the  UTD  expressions  for^  . and  A s !«>  »"d 
(56)  will  be  given  below;  It  wlirhe  shown  that  tncse  expressions  are  valid  within  thd’transitfall  region 
adjacent  to  the  SSB.  Subsequently,  It  will  be  shown  how  these  expressions  automatically  simplify  outside  the 
SSB  transition  region  to  reduce  to  those  obtained  by  Keller  In  his  STB  representation.  Theyt,  h and<A  h In 
(55)  and  (56)  are  [14,151: 


* 

h 


and 


A 


/T-  1 

- -Jm/4 

* , [l-F(XL)1  * P'tL) 

UL 

_ 2/V<  h 

, for  the  lit  region 
(51) 


/STOiWJST 


/f 


e’ij 

k ' (l-F(Xd)j  t P.  (t) 
2/7 1 h 


dn'0l>  e-jkt 


, for  the  shadow  region 
(53) 


TABLE  U 


Zeros  of  the  Airy  Function 

Zeros  of  tneHerivatlve 
of  the  Airy  Function 

A1(-qp)*0 

Ai'(-qp)  * 0 

qj  * 2.33811 

qj  * 1.01879 

q2  - 4,08795 

q,  * 3.24820 

AlM-qj)  * 0.70121 

A1 (-qj)  ■ 0.53566 

Ar(-q2)  - -0.80311 

A1(-q2)  * -0  ‘1  02 

The  function  F appearing  above  has  been  defined  earlier  In  (44),  The  Fock  type  surface  reflection  function 

a soft  P* 

Ps  Is  related  to  the  (hard)  pekeris  function  (q.)  by 
h 
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x 


(Note  that  J*0  at  SSB) 


(59a) 

(59b) 


where  p*  and  q*  are  finite  and  well  behaved  even  when  S*0;  these  universal  functions  are  plotted  In  Figures 
13,  14  and  15.  Also, 


-J; 


\(«'  ■ ~7T~  I dt  e ; 0 . 


MiL  -J5x 


0W2(r) 


r'  r* 

t?  , hard  case 
In  which  the  Fock  type  Airy  functions  V(t)  and  Hj(t)  are 


l , soft  case 


2JV( t ) * w.(t)  - w2(x)  ; W. (r)  • )=  \ dt  eTt't>/3  : 

- e*^,/3 

Vt)  * 75  i j2./3  dt  «rt  1 ^ 

• e 

The  rest  of  the  quantities  occurring  In  (57)  and  (58)  are: 

L , Of  n(t‘)  | kP9(,T|  1/3 

tL  * -ZnlO^Jcose’  ; t ■ / dt  * ^TlfFy  1 n(*)  ■ | ~ — | 


t - /%f  i XL  * 2kL  COsV  ; Xd  • 


(60a) 

(60b) 


(61a;61b) 


(61c) 


(62;63:64) 

(6S;66;67) 


The  quantity  pg  (q^)  In  n(0q)  denotes  the  surface  radius  of  curvature  at  0q  In  the  plane  of  Incidence; 

whereas,  Pg  (Q^)  is  the  surface  radius  of  curvature  at  Qj  In  the  tj  direction.  The  dt'  In  (63)  and  (65)  Is  an 

Incremental  arc  length  along  the  surface  ray  path.  The  angle  of  Incidence  9f  Is  shown  In  Figure  6. 

Also,  the  dn(0j)  and  dn(0?)  In  (58)  denote  the  widths  of  the  surface  ray  tube  at  0i  and  02,  respectively;  the 
surface  ray  tube  Is  formed  by  considering  a pair  of  rays  adjacent  to  the  central  ray  as  In  Figure  11.  The 
geodesic  surface  ray  paths  are  easy  to  find  on  cylinders,  spheres,  and  cones.  For  example,  the  geodesic  paths 
on  a convex  cylinder  are  helical;  whereas,  they  are  great  circle  paths  on  a sphere.  For  more  general  convex 
surfaces,  the  geodesic  surface  ray  paths  must  be  found  numerically.  The  distance  larameter  L In  (66)  and  (67) 
is  given  by 


p}(0,)  pj(0,) 

1 " UjTo^sTTpItVn) 

wheret 


S(PjJ(0j)tS) 

~w~' 


s 


Ph(°l) 


Incident  wavefront  radius  of  curvature 
In  the  bj  direction  at  0j 


(68) 


(69; 70) 


The  distance  s In  (68)  may  be  obtained  by  projecting  (*d)  on  the  SSB  If  the  observation  point  within  the 

(shadow>  stde  of  trans1t1on  region  does  not  rove  In  a predetermined  nanner.  If  the  observation  point 

move'racross  the  SSB  In  a predetermined  fashion  then  It  is  clear  that  s in  (68;69)  , in  he  found  unambiguously. 

The  pj(l)j)  and  p*(0j)  In  (68)  denote  the  principal  radii  of  curvature  of  the  incident  wavefront  at  0j,  and 

pL  which  Is  defined  in  (70),  has  been  introduced  earlier  In  (18b).  For  the  special  case  of  paint  source  or 
spherical  wave  Illumination,  the  l in  (66)  and  (67)  simplifies  to; 


l 


s's 

s‘+s 


where: 


for  spherical  wave  ll'umination, 


s'  = (pjfOj)  - p’(0j)  » p'IOj)) 


distance  from  the  point  source  to  the 
point  of  grating  Incidence  at  0j. 


(71) 

(72) 


In  the  case  of  plane  wave  illumination,  s’ 


and  hence  (71)  above  simplifies  to: 
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L * s , for  plane  wave  Illumination.  (73) 


If  the  Incident  wavefront  Is  of  the  converging  p 
then  the  parameter  L In  (68)  can  become  negative.  It 
solution  can  be  completed  If  L becomes  negative.  On  the  other  hand.  If  one  of  the  principal  directions  of  the 
Incident  wavefront  coincides  with  one  of  the  principal  planes  of  the  surface  at  grating,  then  one  can  treat  a 

converging,  or  converging-diverging  (saddle)  type  wavefront  for  which  L<0,  by  replacing  F(x'"*d)  with 

F*(|X*-,d|),  thjt  the  asterisk  on  F*  denotes  the  complex  conjugate  operator.  The  use  of  F*(|X*‘,d|)  when 

L<0  leads  to  a continuous  total  field  at  SSB  In  this  case. 

The  above  UTD  result  remains  accurate  outside  the  paraxial  (l.e.  near  axial)  regions  of  quasi -cylindrical 
or  elongated  convex  surfaces;  a different  solution  Is  required  In  these  regions  and  It  has  not  yet  been 
completed.  It  Is  assumed  that  the  source  and  observation  points  are  not  too  close  to  the  surface.  Also,  It 
1$  assumed  that  any  caustics  of  the  incident  ray  system  are  not  too  close  to  the  surface.  Furthermore,  the 
amplitude  of  the  Incident  field  Is  assumed  to  be  slowly  varying  at  and  0,;  otherwise.  It  Is  necessary  to 
add  a slope  diffraction  contribution.  The  IITO  solution  described  above  remains  accurate  If  kL  and  m are 
sufficiently  large.  Typically  kL  should  be  larger  than  3 although  In  some  cases  kL  can  be  made  smaller.  Also 
m should  be  such  that  2m3>5  or  so;  however,  the  results  generally  lose  their  accuracy  slowly  as  2m3  becomes 
smaller.  It  Is  noted  that  the  angular  extent  of  the  SSB  transition  region  Is  of  order  m-l  radians. 

A surface  diffracted  field  of  the  type  Ed(P^)  can  also  he  present  In  the  lit  zone  If  the  surface  is 

closed;  this  may  he  seen  by  noting  that  the  field  cf  the  type  Ef*  can  propagate  around  the  closed  surface. 

Also,  additional  contributions  to  (P  ) can  be  present  In  'he  shadow  zone  for  a closed  surface  because 
surface  diffracted  rays  can  be  Initiated  at  all  points  of  grazing  Incidence  on  that  closed  surface; 
furthermore,  these  surface  rays  can  undergo  multiple  encirclements  around  the  closed  body.  However,  these 

additional  surface  diffracted  ray  contributions  are  generally  quite  weak  in  comparison  to  the  F^r  contribution 
within  the  lit  zone  for  surfaces  which  are  quite  large  In  terms  of  the  wavelength;  hence  their  contribution 
may  be  neglected  In  such  cases. 

The  parameters  t*",  t,  X*-  and  x'1  become  small  as  one  approaches  the  surface  shadow  boundary,  SSB,  from 
both  the  lit  and  shadow  regions.  As  one  approaches  the  SSB,  the  small  argument  limiting  form  of  the 
transition  function  F(X)  which  has  been  introduced  previously  In  (44)  becomes  helpful  for  verifying  the 
continuity  of  the  total  high  frequency  field  at  the  SSB.  On  the  other  hand,  the  above  parameters  become  large 
as  one  moves  outside  the  SSB  transition  region;  In  this  case  F(X)*1  for  large  X,  and  likewise, 


! - < 0,  or  convorglng-dlverglng  (p!  < 0;  pi  > 0)  type, 
has  not  been  fully  Investigated  how  the  general 


?,(«> 

h 


a«o 


Pj(4) 

h 


4>>0 


( 


■J5./6 


n l 2[A1 1 !-q  ))2 ” 
n*l  n 


-j»/4 


Ji/6  Jq  e 


■J5./6 


e ' N e e 

7 r*  i 


n»i  «VMK>V 


(74) 


(75) 


where  N * 2 is  generally  sufficient  to  compute  ^5(4)  accurately  for  4»0  In  (75).  In  (7b)  and  Table  II,  the 

Hiller  type  Airy  function  A1(t)  * V(t)/V»,  and  A 1 ‘ ( t ) -g-AI(t),  Thus,  upon  incorporating  the  limiting 
values  of  f 74 } and  (75),  which  are  valid  outside  fhe  SSB  jransltlo  • c.jion,  Into  (57)  and  (58)  and  replacing 
F(X)  by  It  * s asymptotic  value  of  unity,  ft  Is  clear  (hat7ls  revues  to  <?s  * ? 1 outside  the  SSB  transition 

region  so  that  E^r (P^ ) < Er(P^)  of  GO,  and  likewise  Ed(Ps ) > £b(Ps)  therein,  respectively.  In  which  the  Keller 

surfaced  diffracted  ray  field  Ed(Psl  is  given  by  [3) 


E?(s) 


e'(o,)- 


k(0j,02)  e 


■jkt 


/dn(Qi)  F 
/ dn(Qj)  / 


<>s 


sd(p5nd> 


-jkt 


(76) 


Tk<W  * I h\  \ Ts  * "1  "2  Tn 


where 


TS  * i- 


02  s 

-/  « " (t')df 

Ql  " 


n=»l 


o*  (Op) 


(77) 


tOi ) e 


(78) 
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The  ojj  and  « jj  are  the  Keller's  GTO  diffraction  coefficients  and  attenuation  constants  for  the  nth  soft  (s)  or 

hard  (h)  surface  ray  mode.  Thus,  In  the  GTO,  the  surface  ray  field  consists  of  surface  ray  modes  which 
propagate  Independently  of  one  another.  Also,  this  surface  ray  field  Is  not  the  true  field  on  the  surface;  It 

Is  a boundary  layer  field.  The  0^  and  « 8 are  given  by: 

n n 

.-j(*/12) 


[0„${Q)]2  * /£ 


"(Q)  [A1'(-q„)]t 


(79a) 


i 


h 2 rr  * 
[0„(Q)1  -/m  «(q)  — 


-J<»/12) 


9nW1  (-Rn^ 


and 


s #A,  In  , h 5n  , J(*/6) 

« n W)  " pJWT  ■«>  : * n 


(79b) 


(80a;80b) 


In  (79)  - (80),  <)  Is  any  point  on  the  geodesic  surface  ray  path.  The  GTD  result  of  (76)  In  terms  of  (77)  and 
(78)  Is  not  valid  within  the  SSB  transition  region. 

The  UTD  result  for  the  3-0  configuration  can  be  simply  modified  to  recover  the  corresponding  UTO  result 
for  the  2-0  case  by  allowing  the  caustic  distances  nj  and  o5  In  (55)  and  (56)  to  receed  to  Infinity.  Then, 


P P P 

pj  i s , If  uj  ♦ « and  p$ 


so  that 


csr(pL)  ^ e1^).  {fts  ;A  ;A  ♦ x*  ;j, //  e-Jk*r 

V pp+sr 

In  which  pr  Is  as  In  (23),  and 

fd(%)  - E1  (0,).  (A  b b + A n,  n?}  e ^ks 

TP 


(81) 


(82) 


(83) 


■R  s and  A 
h 

respectively;  only  the  L appearing  In  (66)  and  (67)  Is  given  by 


since  (bj  » b?  z b for  the  2-0  case  (note:  b » e^)).  The  As  and  A In  (82)  and  (83)  are  as  defined  earlier, 

h h 


s'sd 


l ■ s.+sd  , for  the  2-0  case, 


(84) 


i where  s'  Is  the  distance  from  the  2-11  line  source  to  the  point  of  grazing  Incidence  at  0,  and  s = sd  as 

1 SSB 

l before.  A comparison  of  the  UTO  and  GTO  solutions  for  a 2-0  circular  cylinder  Illuminated  by  a nearby  line 

! source  is  Illustrated  in  Figures  16(a)  and  16(b);  those  UTD  solutions  are  then  compared  with  the 

j corresponding  exact  (Eigenfunction)  solutions  in  Figures  17(a)  and  17(b). 

III.  A FEU  ADDITIONAL  UTD  SOLUTIONS  AND  SOME  APPLICATIONS 

| 

In  addition  to  the  basic  UTD  edge  and  convex  surface  diffraction  solutions  described  above,  UTD  solutions 
for  some  other  canonical  shapes  also  exist;  however,  the  latter  are  not  described  here  because  of  space 
limitations.  UTO  type  solutions  for  the  radiation  and  mutual  coupling  associated  with  antennas  on  a smooth 
convex  surface  are  given  In  T16-19];  also,  an  approximate  vertex  diffraction  solution  may  be  found  In 
[5,11,203.  A result  based  on  a recently  obtained  approximate  UTO  solution  for  the  field  scattered  by  a fully 
Illuminated,  semi-infinite,  right -circular  perfectly-conducting  cone  [?H  is  shown  In  Figure  18.  Also,  UTO 
results  for  the  3-D  diffraction  by  a penetrable  dielectric/ferrite  s.rip  in  Figure  19  based  on  the  work  In 
[221  are  shown  In  Figures  20  and  21,  for  parallel  and  perpendicular  polarization  of  the  Incident  field, 
respectively.  It  Is  noted  therein  that  even  though  the  Incident  fields  are  TEZ  or  TXZ,  the  scattered  fields 
are  not  simply  TEZ  or  THj  due  to  a coupling  between  the  two  which  is  introduced  by  the  dielectric  edge  when 
®'*»/2.  Finally,  Figures  22  and  23  show  the  application  of  UTO  to  deal  with  more  realistic  shapes  [23,241. 

The  ogival  shape  In  Figure  22  has  a circular  duct  on  It.  In  Figure  23,  the  aircraft  fuselage  Is  modeled  by  a 

best  fit  prolate  spheroid  near  the  antenna  location,  a mare  recent  calculation  employs  a composite  e'lfpsoid  % 

fuselage  model  [25]. 
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(b)  diffraction  at  a discontin- 
uity In  surface  Impedance 
(Zsx  and  ZS2> 


(c)  diffraction  at  a discontin- 
uity In  surface  curvature 


(d)  diffraction  by  a thin  curved 
wire 


(e)  diffraction  by  a smooth 
convex  surface 


(f)  diffraction  by  a vertex  In  a 
plane  screen 


Figure  t.  Examples  of  diffracted  rays 
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Figure  2,  Equlphase  (or  wavefront)  surface 


Figure  5.  Reflected  wavefront  geometry 
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Figure  3.  Ray  tube  geometry  associated  with 
the  ray  path  TSP 
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Figure  4.  Examples  of  shadow  and  reflection 
boundaries 


Figure  6.  Unit  vectors  fixed  In  the  Incident 
and  reflected  rays 


Figure  7.  Edge  diffracted  ray  tube 
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16.  Comparison  between  GTO  and  UTD  solutions  for  the  radiation  by  electric  (a)  and  magnetic  (b)  line 
sources  In  the  presence  of  a circular  cylinder.  Here  a*lx  and  p'»2X. 
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Figure  17,  Comparison  of  UTO  solutions  of  Fig,  16  with  exact  (eigenfunction)  solutions 


18,  Magnetic  field  Ha  scattered  by  a semi-infinite,  perfectly-conducting  right  circular  cone  when 
illuminated  with  a plane  wave  e! 


